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The dynamical dielectric function of two dimensional graphene at arbitrary wave vector q and 
frequency to, e(q,ui), is calculated in the self-consistent field approximation. The results are used to 
find the dispersion of the plasmon mode and the electrostatic screening of the Coulomb interaction 
in 2D graphene layer within the random phase approximation. At long wavelengths (q — > 0) the 
plasmon dispersion shows the local classical behavior uj c i = uja^/q, but the density dependence of the 
plasma frequency (wo oc 7i 1//4 ) is different from the usual 2D electron system (wo oc n 1 ^ 2 ). The wave 
vector dependent plasmon dispersion and the static screening function show very different behavior 
than the usual 2D case. We show that the intrinsic interband contributions to static graphene 
screening can be effectively absorbed in a background dielectric constant. 

PACS numbers: 71.10.-w; 73.21.-b; 73.43.Lp 



INTRODUCTION 

There has been a great deal of recent interest in the 
electronic properties of two dimensional (2D) graphene, 
a single-layer graphite sheet, both theoretically and ex- 
perimentally [l|, [2| . The main difference of 2D graphene 
compared with other (mostly semiconductor-based) 2D 
materials is the electronic energy dispersion. In conven- 
tional 2D systems, the electron energy with an effective 
mass m* depends quadratically on the momentum, but 
in graphene, the dispersions of electron and hole bands 
are linear near K, K' points of the Brillouin zone. Be- 
cause of the different energy band dispersion, screening 
properties in graphene exhibit significantly different be- 
havior from the conventional 2D systems [3J] . The screen- 
ing of Coulomb interaction induced by many body effects 
is one of the most important fundamental quantities for 
understanding many physical properties. For example, 
dynamical screening determines the elementary excita- 
tion spectra and the collective modes of the electron sys- 
tem, and static screening determines transport properties 
through screened Coulomb carrier scattering by charged 
impurities. In this paper, we theoretically obtain the (dy- 
namical and static) screening behavior of 2D graphene 
by calculating, for the first time, the polarizability and 
the dielectric function within the self-consistent field ap- 
proximation (i.e. random-phasc-approximation (RPA)) 
for gated-2D graphene free carrier systems. We apply 
our theory to calculate the 2D graphene plasmon disper- 
sion and the static screening function, finding some inter- 
esting qualitative differences between graphene and the 
extensively studied 2D electron systems based on semi- 
conductor heterostructures and MOSFETs. 

In this paper, we calculate the dielectric function of 
graphene at arbitrary wave vector q and frequency u>, 
e(q,iv), within RPA, in which each electron is assumed 
to move in the self-consistent field arising from the ex- 
ternal field plus the induced field of all electrons. This 



is the model which leads to the famous Lindhard dielec- 
tric function for a three-dimensional (3D) 0] and 2D Q 
electron gas. One of the immediate theoretical conse- 
quences of the dielectric function is that its zeros give 
the wave vector dependent plasmon mode, u> p i(q), which 
is a fundamental elementary excitation and a collective 
density oscillation mode. Using the theoretical dielectric 
function we provide the plasmon mode dispersion both 
for single-layer and bilayer graphene. Another important 
consequence of the dielectric function is the static screen- 
ing function which can be obtained as the static limit 
ijj — > of the dielectric function, describing the electro- 
static screening of the electron-electron, electron-lattice, 
and electron-impurity interactions. 



POLARIZABILITY: 11+ AND fl 

The electron dynamics in 2D graphene is modeled by 
a chiral Dirac equation, which describes a linear relation 
between energy and momentum. The corresponding ki- 
netic energy of graphene for 2D wave vector k is given 
by (we use h = 1 throughout this paper) 



e sk = S7|k|, 



(1) 



where s = ±1 indicate the conduction (+1) and va- 
lence (—1) bands, respectively, and 7 is a band param- 
eter (essentially the 2D Fermi velocity, which is a con- 
stant for graphene instead of being density dependent). 
The corresponding density of states (DOS) is given by 
D(e) = g s g v \e\/(2iry 2 ), where g s = 2, g v = 2 are the spin 
and valley degeneracies, respectively. The Fermi momen- 
tum (hp) and the Fermi energy (Ep) of 2D graphene are 
given by kp = (47m/ gsgv) 1 ^ 2 and Ep = ^ykp where n is 
the 2D carrier (electron or hole) density. For the sake 
of completeness, we also mention that the dimensionless 
Wigner-Seitz radius (r s ), which measures the ratio of the 
potential to the kinetic energy in an interacting quantum 
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Coulomb system [J], is given in doped 2D graphene by 
r s = (e 2 /K7)(4/g s <7„) 1 / 2 where k is the background lat- 
tice dielectric constant of the system. We note in the 
passing the curious fact that the dimensionless r s pa- 
rameter is a constant in graphene unlike in the usual 2D 
(r s ~ n -1 / 2 ) and 3D (r s ~ ?i -1 / 3 ) electron liquids, where 
r s (and therefore interaction effects) increases with de- 
creasing carrier density. The constancy of r s in graphene 
arises trivially from the relativistic Dirac-like nature of 
the free carrier graphene dynamics implying that the 'rel- 
ativistic' effective mass, m c = Ep/j 2 , depends on carrier 
density precisely as y/n cancelling out the corresponding 
yfn term in the potential energy. Equivalently, r s here is 
just the "effective fine structure constant" for graphene, 
with a value of r s ~ 0.5 assuming g s = g v = 2 and 
k = 4 (using Si02 as the substrate material). This small 
(and constant) value of graphene r s indicates it to be a 
weakly interacting system for all carrier densities, mak- 
ing RPA an excellent approximation in graphene since 
RPA is asymptotically exact in the r s < 1 limit. 

In the RPA, the dynamical screening function (dielec- 
tric function) becomes 



e(q,u>) = l + v c (q)U(q,u), 



(2) 



where v c (q) = 27re 2 / ' nq is the 2D Coulomb interaction, 
and n(q,aj), the 2D polarizability, is given by the bare 
bubble diagram 
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(3) 

where k' = k + q, s,s' = ±1 denote the band in- 
dices and F ss /(k, k') is the overlap of states and given 
by F ss i(\s., k') = (1 + ss' cos 9)/2, where 9 is the angle 
between k and k', and f s k is the Fermi distribution func- 
tion, / sk = [exp{/3(e sk -^)} + l]- 1 , with /3 = l/k B T and 
\x the chemical potential. After performing the summa- 
tion over ss' we can rewrite the polarizability as 



where 
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n-( g ,w) 



g s gv 

~ 2L 2 



E 



[/k--/k'-](l + cos^ fe 



/k-(l - cos 9 k k>) 



lu + e k _ - e k /_ + IT] 
/k'-(l - cos 9 kk ') 



lu + e k _ 



£k' 



IT] 



tu + e kH 



£k'- 



IT} 



(0) 



For intrinsic (i.e. undoped or ungated with n and Ep 
both being zero) graphene, in which the conduction band 



is empty and the valence band is fully occupied at zero 
temperature (i.e. Ep = 0), we have / k + = and / k _ = 
1. Then the polarizability becomes U(q,Lu) = U^(q,Lu), 
which has been previously obtained in the renormaliza- 
tion group approach Q. Recently Tl~(q,Lu) has been 
reconsidered to discuss screening effects of Coulomb in- 
teraction in intrinsic graphene [7[. In general, H + (q,Lu) 
does not vanish for most systems because the Fermi en- 
ergy is typically located in the conduction or the valence 
band. But graphene is a most peculiar zero-gap semi- 
conductor system where Ep — in the intrinsic undoped 
situation. In the doped or gated situation n, Ep ^ 
in graphene, and now II + is finite. In the following we 
provide the zero temperature polarizability in the doped 
or gated case where the Fermi energy is not zero. 

By introducing the dimensionless quantities x = q/kp 
and v = Lu/Ep, and tl([q,u>) = H([q,Lu)/ Dq where Dq = 
D(E F ) = (gsgyn/n) 1 / 2 /"/ is the DOS at Fermi energy, 
we have 

n+(:r, v) = fLf(x, v)9(y - x) + ft+(x, v)9{x - v), (7) 

where the real parts of the polarizability are given 

1 
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+ —[e(\v + 2\-x) + e(\v-2\-x)}}(9) 
and the imaginary parts of the polarizability are 
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where 
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f 3 (x, v) = (2 + v) v/x 2 - (2 + vf + x 2 sin- 1 (14) 
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and II (x,f) can be calculated to be 



n (x,v) = — , i J + 1 — , I ' (16) 
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Eqs. (|7)l- (fTB|) are the basic results obtained in this paper, 
giving the 2D doped graphene polarizability analytically. 
Note that our 2D graphene polarizability is completely 
different from the corresponding 2D Lindhard function 
first calculated in ref. la, which is appropriate for the 
usual 2D systems with parabolic band dispersion. 



PLASMONS IN RPA 

As a significant consequence of the dielectric func- 
tion we calculate the long wavelength plasmon disper- 
sion for single-layer graphene and for bilayer graphene. 
The longitudinal collective-mode dispersion, or plasmon 
mode dispersion, can be calculated by looking for poles 
of the density correlation function, or equivalently, by 
looking for zeros of the dynamical dielectric function, 
e(q, uj) = 1 — w(g)IT(q, uj). In the long wavelength limit 
(q — > 0) we have the following limiting forms in the high- 
and low-frequency regimes: 
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, 7g < uj < 2Ep 
uj < jq. 



(17) 



In the q — > limit, we have the plasmon mode dispersion 
uj p (q) for a single-layer graphene as 



uj c i = uj p {q -> 0) = ujo^/q 



(18) 



where ujq = [g s g v e 2 Ep 12k) 1 ! 2 . The leading order (or 
local) plasmon has exactlythe same dispersion, q 1 ^ 2 , as 
the normal 2D plasmon [3j. However, the density de- 
pendence of the plasma frequency in graphene shows a 
different behavior, i.e. ujq oc n 1 / 4 compared with the 
classical 2D plasmon behavior where ujq oc n 1 / 2 . This is 
a direct consequence of the quantum relativistic nature 
of graphene. Even though the long wavelength plasmons 
have identical dispersions for both cases, the dispersion 
calculated within RPA including finite-wave-vector non- 
local, (i.e. higher order in q) effects show very different 
behavior. In normal 2D 0, 0| the non-local correction 
leads to an increase in plasma frequency, [u> p (q)/uj c i = 
1 + (3/4)(q/qTF)], where qxp = Ssff-u^e 2 //* is the usual 



2D Thomas-Fermi wave vector) , but in graphene the cor- 
rection within RPA leads to a decrease in plasma fre- 
quency compared with ui c i [uj p (q)/uj c i = 1 — q q/8kp], 
where go = 9s9v^ 2 ^f ll^ is the corresponding graphene 
Thomas-Fermi wave vector. Recently, the plasmon mode 
of graphene has been considered numerically in the pres- 
ence of spin-orbit coupling Q. 

For bilayer graphene, without inter-layer hopping, we 
have the leading order q dependence of the collective 
modes by solving a two component determinantal equa- 
tion 



uj-(q) 



uj oy J2q 
2uj Vdq, 



(19) 



where d is the layer separation between the two 2D 
graphene sheets. The w+ mode, the optical-plasmon 
mode (in-phase mode of the coupled system), has the well 
known q 1 / 2 behavior, independent of the layer separation 
d at long wavelengths. The other mode uj_ is the acoustic 
plasmon mode (out-of-phasc mode of the coupled system) 
which goes as q in long wavelengths and depends on the 
separation d. Thus, the coupled plasmons in graphene 
show the same long wavelength behaviors as those of 
normal 2D systems. But, the density dependences of 
the plasma frequency and the large wave vector depen- 
dences are again very different from the corresponding 
normal 2D systems [9j. When interlayer hopping is in- 
cluded in a bilayer system, the in-phase plasmon mode 
is qualitatively unaffected by tunneling. However, the 
out-of-phasc plasmon mode develops a long wavelength 
gap (depolarization shift) in the presence of tunneling 
[TJ], i.e. ui-(q — > 0) = (2t±) 2 (l + qod), where t± is the 
inter-layer hopping. Due to strong interlayer coupling in 
bilayer graphene, in general, we have uj+ <^ U- at long 
wavelengths. 

In Fig. Q] we show the calculated plasmon dispersion 
within RPA (solid line) compared with the classical local 
plasmon (dashed line). We use the following parameters: 
k = 2.5, 7 = 6.5eVA, and a density n = 10 12 cm -2 . 
In fig. [TJb) we show the corresponding 2D regular plas- 
mons with n-GaAs parameters (n = 5 x 10 11 cm' 2 ). In 
Fig. [T]we also show the electron- hole continuum or sin- 
gle particle excitation (SPE) region in (q, uj) space, which 
determines the absorption (Landau damping) of the ex- 
ternal field at given frequency and wave vector. The SPE 
continuum is defined by the non-zero value of the imag- 
inary part of the polarizability function, lmll(q,uj) ^ 0. 
For a normal 2D system only indirect (q ^ 0) transi- 
tion is possible within the band, and the SPE bound- 
aries are given by u>i t 2 = q 2 /2m ± qkp/m. However, for 
2D graphene both intraband and interband transition are 
possible, and the boundaries are given in Fig. 1(a). The 
intraband SPE boundaries are uj\ = jq (upper bound- 
ary) and UJ2 — for q < 2kp, uii = 7? — 2Ep for q > 2Uf 
(lower boundary). The direct transition (q = 0) is also 
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possible from the valence band to the empty conduction 
band. Due to the phase-space restriction the interband 
SPE continuum has a gap at small wave vectors. For 
q = 0, the transition is not allowed at < u> < 2E F . 
If the collective mode lies inside the SPE continuum we 
expect the mode to be damped. Since the normal 2D 
plasmon lies, at long wave lengths, above the SPE con- 
tinuum it never decays to electron-hole pair within RPA. 
But for graphene the plasmon lies inside the interband 
SPE continuum decaying into electron-hole pairs. Only 
in the region I of Fig. 1 (a) the plasmon is not damped. 
The other different feature between a normal 2D plas- 
mon and a graphene plasmon occurs at large wave vec- 
tors. The normal 2D plasmon mode enters into the SPE 
continuum at a critical wave vector, and therefore does 
not exist at very high wave vectors. All spectral weight 
of the plasmon mode is transferred to the SPE. But the 
graphene plasmon does not enter into the intraband SPE 
and exists for all wave vectors, except for its decay into 
real interband electron-hole pairs in the SPEi ntor regime. 

STATIC SCREENING 

Now we consider the static polarizability H(q, u> = 0). 
From Eq. (O we have 



fi+(g) 



1 
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and from Eq. (|16[) we have 



q < 2k F 

r 1 ^, q>2k F 
(20) 



fl~{q) = nq/8k F . 



Thus, the total static polarizability becomes a constant 
at q < 2k F as in a normal 2D systems, i.e. H(q) = 
IL+(q) + U-(q) = D(E F ) for q < 2k F . In Fig. |wc 
show the calculated static polarizability as a function of 
wave vector. For a normal 2D system the screening wave 
vector, q s = q FF = g s g v me 2 / 'k, is independent of elec- 
tron concentration, but for 2D graphene the screening 
wave vector is given by q s = g s g v e 2 k F / kj which is pro- 
portional to the square root of the density, n 1 / 2 . In the 
large momentum transfer regime, q > 2k F , the static 
screening increases linearly with q due to the interband 
transition. This is a very different behavior from a nor- 
mal 2D system where the static polarizability falls off 
rapidly for q > 2k F with a cusp at q — 2k F [§]. The 
linear increase of the static polarizability with q gives 
rise to an enhancement of the effective dielectric constant 
n*(q — > oo) = k(1+ g s g v Trr s /8) in graphene. Note that in 
a normal 2D system k* — ► k as q — ► oo. Thus, the effec- 
tive interaction in 2D graphene decreases at short wave 
lengths due to polarization effects. This large wave vector 
screening behavior is typical of an insulator. Thus, 2D 
graphene screening is a combination of "metallic" screen- 
ing (due to n + ) and "insulating" screening (due to II - ), 



leading to overall rather strange screening properties, all 
of which can be traced back to the zero-gap chiral rela- 
tivistic nature of graphene. 

In may be worthwhile to ask whether the intrinsic 
graphene contribution, arising strictly from the interband 
transitions due to the filled valence band (i.e. the II - 
term in our graphene polarizability), can be absorbed 
in the effective background lattice dielectric constant K 
just as one does in a regular semiconductor (usi = 11.5; 
KGaAs = 12.9) or insulator (nsi0 2 = 3.9) in discussing 
free carrier screening by doping or gating indeed free car- 
riers in conduction (electrons) or valence (holes) bands. 
In particular, only intraband free carrier screening is ex- 
plicitly considered in the usual 2D screening function 
[H, HI extensively used HI, 3] in the quantitative analy- 
sis of quantum transport in 2D semiconductor devices 
such as Si MOSFETs, GaAs modulation-doped high- 
mobility transistors (HEMTs), and undoped gated GaAs 
heterostructures (HIGFETs). The interband transition 
induced screening in the semiconductor based 2D struc- 
tures is included in the theory simply by appropriately 
modifying the effective background lattice dielectric con- 
stant from the usual vacuum value of unity to a value 
around ten. 

To see whether the effect of the interband n polariz- 
ability can be 'trivially' absorbed in a background lattice 
dielectric constant we rewrite a 2D graphene dielectric 
function e(q) = 1 + v c (q)H(q) to obtain: 



2ttp 2 

s(q) = l + [n-(q)+U+(q)} 

Kq 



(22) 



(21) Using Eq. (EU), U~(q) = D{E F )Trq/8k F} we get 



« c (?)H-(g) = ^f^^ = ^f^r s . Thus, we have 

e(q) = l + i^Lr s +v c (q)U+(q). (23) 

Introducing an effective intrinsic background graphene 
dielectric constant 



we have 



e{q) = n* 



1 + g s g v -KT s 



2ne 2 , , / 
KK*q 



(24) 



(25) 



Writing an effective free carrier 2D graphene dielectric 
function 



e + (q) = l+v+(q)U+(q), 
where v^(q) = 2ire 2 / KK*q, we have 
e(q)^K*e + (q). 



(26) 



(27) 



Eq. |27|) shows that the intrinsic screening contribu- 
tion arising from the interband n~ term can be com- 
pletely subsumed by introducing an effective graphene 
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background lattice dielectric constant k* = l-\-g s g v 7rr s /8, 
and by making the replacement k — ► kk* throughout. In- 
troducing the effective dielectric constant k* allows one 
to use only the free carrier screening function 



1 



KK*q 



n + (<z), 



(28) 



for describing free carrier screening properties of 2D 
graphenc. We note that k itself here is the back- 
ground lattice dielectric constant arising from the in- 
sulating substrate (i.e. SiC>2 in most situations) with 
k = (1 + Ksio 2 )/2 ~ 2.5, and k* = 1 + g s g v irr s /8 fw 2.3 



with 



e / nh/y « 0.7 for graphene on Si02 substrate. 



Thus, the substitution k — * kk* arising from the in- 
tcrband contributions enhances the effective background 
graphenc dielectric constant to kk* sa 6 — this approxi- 
mate factor of 2 increase in the background dielectric con- 
stant arising from interband contributions further sup- 
press Coulomb interaction effects in extrinsic graphenc. 

The large q behavior of graphene dielectric screening, 
e{q — > oo) — > k* = 1 + g s g v irr s /8, again demonstrates 
that the short-wavelength screened Coulomb potential in 
2D graphene goes as 2ire 2 / KK*q, with an enhanced back- 
ground effective dielectric constant kk* where k is the 
effective background dielectric constant arising from the 
substrate and k* arising from graphene interband po- 
larizability II as discussed above. This implies that a 
suspended 2D graphene film, without any substrate (i.e. 
k = 1), would have an effective background lattice dielec- 
tric constant of k* = l+g s g v irr s /8 with r s = e 2 /fvy (since 
k = 1). Putting in g s g v = 4, we get k* rs 4. Thus the 
background static lattice dielectric constant of intrinsic 
graphene, due to interband transitions, is around 4. 

If the T = transport properties of graphenc arc 
dominated by c harg ed impurity scattering as is thought 
to be the case 15j then the long-wavelength Thomas- 
Fermi screening becomes an important input for calcu- 
lating the screened charged impurity potential: €tf{q) = 
£rpa(q — ► 0) becomes 



£tf(<?) = 1 + <7tf/<?, 



(29) 



where qTF = <Z« = gsgvG 2 kp j kj. Note that one can 
equivalently define a long wavelength effective TF screen- 
ing function e^, p {q) where interband screening effects are 
absorbed in an effective background dielectric constant 



£tf(i) = 1 + <7tf/<7, 



(30) 



where q^ F = qTF / K * ■ As discussed above, the identity, 
e(q) = K*e + (q), guarantees the equivalence between Eqs. 
(|2T)1) and (|3"0"|) for long wavelength graphenc screening. 
We note also that the q = screening wavevector qTF is 
simply proportional to the graphenc density of states at 
the Fermi level at T = 0. 



Although the two screening descriptions, based on e(q) 
with the background dielectric constant being just k and 
on e + (q) with the background dielectric constant being 
kk* , are precisely equivalent for T = static screening 
properties, the two descriptions are not inequivalent at 
finite temperatures. Therefore, it is more appropriate to 
use the full dielectric function e(g) in theoretical work on 
2D graphenc. 



CONCLUSION 

In conclusion, we have theoretically obtained ana- 
lytic expressions for doped (i.e. Ep ^ 0) 2D extrinsic 
graphenc polarizability, dielectric function, plasmon dis- 
persion, and static screening properties, finding a num- 
ber of intriguing qualitative differences with the corre- 
sponding normal (and extensively studied) 2D electron 
systems. The differences, with interesting observable 
consequences, can all be understood as arising from the 
zero-band gap intrinsic nature of undoped graphene with 
chiral linear relativistic bare carrier energy band disper- 
sion. Some of our qualitatively new predictions, such as 
the n 1 / 4 dependence of the long-wave length graphenc 
plasma frequency in contrast to the well-known n 1 / 2 be- 
havior of classical and normal 2D plasmons, should be 
easily verifiable experimentally using the standard exper- 
imental techniques of infra-red absorption and/or in- 
elastic light scattering [T2| spectroscopies. Similarly, our 
prediction of the peculiar nature of the graphenc plas- 
mon damping (i.e. no Landau damping due to intraband 
electron- hole pairs, but finite Landau damping due to in- 
terband electron-hole pairs) should be easily verifiable. 
Our predicted different screening behavior in graphenc 
at large wave vector should have consequences for trans- 
port properties. Our RPA theory should be an excellent 
qualitative approximation for 2D graphenc properties at 
all carrier densities (as long as the system remains a ho- 
mogeneous 2D carrier system, which may not be true 
for n < 10 12 cm~ 2 ) since the effective rvparameter for 
graphene is a constant (< 1), making RPA quantitatively 
accurate in graphene. Finally, we point out that the ef- 
fective Fermi temperature, Tp = Ep/ftB, being very high 



(~ 1300AT for 



10 cm ) in graphene, our T = 



theory should apply all the way to room temperatures. 
We note that the long wavelength dielectric function for 
bulk graphite was earlier considered within an approxi- 
mation scheme in ref. [Til and the zero frequency limit 
was recently considered in ref. [17. 

Before concluding, we point out that some aspects of 
graphene collective modes and linear response have been 
discussed in the recent literature. In particular, the in- 
trinsic situation without any free carriers has been con- 
sidered in [3] whereas our emphasis in this work has been 
extrinsic graphene with free carriers (electrons/holes in 
conduction/ valence band) induced by external gating or 
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doping. There has been a recent purely numerical study 
[1] of graphene collective mode spectra in the presence of 
spin-orbit coupling. 

Note added After submitting the present paper, we be- 
came aware of related work [18( . 

This work is supported by the US-ONR, the LPS, and 
the Microsoft corporation. 
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FIG. 1: (a) Plasmon mode dispersion in 2D graphene (solid 
thick line) calculated within RPA. The dashed line indicates 
the local long wavelength plasmon dispersion. Thin solid 
lines represent the boundaries of the single particle excita- 
tion (SPE) Landau damping regime for intra- and inter-band 
electron-hole excitations, (b) The plasmon dispersion and 
SPE for a normal 2D system with a quadratic energy disper- 
sion. 




FIG. 2: Static polarizability for 2D graphene: Utot = n + 
II". Here 11(0) = D(E F ) = g s g v k F /2tvj. 



